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This study investigates the feasibility of applying the solution 
techniques of Mlnlquasi to the set of equations which govern the 
thermoregulatory model of Stolwijk. For solving nonlinear equa- 
tions and/or boundary conditions, a Taylor Series expansion is 
required for linearization of both equations and boundary con- 
ditions. The solutions are iterative and in each iteration, a 
problem like the linear case is solved. It will be shorn that 
Miniquasi cannot be applied to the thermoregulatory model of 
Stolwijk as originally planned. 
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PROGRAM DESCRIPTION GUIDE 


A. IDENTIFICATION 

Program Name - MINIQUASI and the Thermoregulatory Model of 

Stolwijk 

Programmer's Name - V. J. Marks 
Date of Issue - September 20, 1974 


B. GENERAL DESCRIPTION 

This study will investigate the feasibility of applying the solution 
techniques of Miniquasi to the set of equations which govern the ther- 
moregulatory model of Stolwijk. For solving nonlinear equations and/or 
boundary conditions, a Taylor Series expansion is required for lineari- 
zation of both equations and boundary conditions . The solutions are 
iterative and in each iteration, a problem like the linear case is 
solved. It will be shown that Miniquasi cannot be applied to the thermo- 
regulatory model of Stolwijk as originally planned. 


C. USAGE AND RESTRICTIONS 

Machine and Compiler Required 


Peripheral Equipment Required 

Approximate Amount of Memory 
Required 


D. PARTICULAR DESCRIPTION 

Equations Used and Derivations - See References 1 and 2 

Definitions of Terms Used - See Reference Set of 1 - Appendix A 

and B 

The Miniquasi program is designed to provide numerical solution of- multi- 
point boundary value problems governed by a system of up to 20 first- 
order ordinary differential equations. This system of equations may have 
up to 20 boundary conditions . The equations must be written .and pro- 
grammed in their linearized form for purposes of applying Miniquasi. The 
equations which describe the controlled system of the thermoregulatory 
model are given below with temperature as the dependent variable T . 

"fctl o 

T(i) Temperature of the l element, F. 

t Time, Hours 

VPDEW Vapor Pressure at Cabin Temperature, mm Hg 


- Xerox Sigma 3 

- ANSI FORTRAN 

- Card Reader and Line Printer 

- 1382 hexidecimal (504-3 decimal) 
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C(x) 

Specific Heat * Mass, BTU/°F 

P 

Cabin Pressure, PSI 

TCAB 

Cabin Temperature, °F 

TUG(i) 

Temperature of Garments, °F 

TW 

Temperature of the Cabin Wall, °F 

CLO 

Ratio of Clothed Area to Total Area 

RM 

Metabolic Rate, BTU/Hr.' 

QB(i) 

Basal Metabolic Rate, BTU/Hr. 

UEFF 

Useful Work Efficiency, $ 

AC 

2 

Total Convective Area, Ft 

AR 

2 

Total Radiative Area, Ft. 

VCAB 

Cabin Free Stream Velocity, CFM 

VEFF 

Cabin Ventilation Efficiency, °lo 

G 

Gravity (Normal to Earth) 

CLOV 

Clothing Thermal Resistance, Hr. Ft 

EUG 

Emissivity of Undergarment 

CPG 

Specific Heat of Gas, BTU/lb. °F 


°f/btu 


The thermoregulatory model contains 4l nodes composed of the core, muscle, 
fat, and skin for the ten geometric compartments of the head, trunk, arms, 
legs, hands, and feet. The forty-first node is the central blood pool. 

The equations of the heat balance in these compartments are of the basic 
form. 

Mass * Cp * H =AQ (!) 

Where A Q, represents the heat loss or heat production terms associated 
with each of the forty-one nodes . The terms include conduction to adja- 
cent nodes, convection to the blood, radiation terms, and the sensible and 
latent heats of vaporization from both the skin and the lungs . The many 
assumptions necessary to derive the heat flow terms from physical laws and 
geometry can be found in References 2 and 3- 
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It is the heat terms when expanded into functions of only independent 
variables which are nonlinear and are not amenable to the linearization 
techniques suggested in Miniquasi. These equations in their expanded 
form become: 


Temperature of the Head Core = T(l) 


[• 


dT(l) = 1 * QMET(l) + Q,C0NV(l) - QELATl 

at e(i) 


( 2 ) 


-Q,C0ND(l) - QpSEN(l)J 


where these terms are from 

QMET(i) Metabolic heat produced in the i™ 1 element. 

QCONV(i) Heat exchanged due to convection between blood and 

tissue in the i^“ element . 

QRLAT1 The component of latent heat loss from the lungs which 

is distributed to the head core . 

QCOUD(i) The heat exchanged with adjacent tissue compartments 

due to conduction in the x element . 

QRSEN(i) The sensible heat losses from the lungs which are dis- 

tributed to the l element. 


In expanded form: 

QMET(l) = QB(l) =( CONSTANT) (3) 

QCONV(l) = BF(l) * (T(4l) - T(l) ) (4) 

where 

BF(l) is the blood flow rate in the compartment ? &nd 

BF(l) = 99'. 3 lb/hr. (5) 

0.386 * 0 . 0415 * PCAB * 144 .0 
QKLAT1 - 51.5 * (TCAB + 460.) 

* [l.O - .33 *(14.7 - PCAB)| * EM 

* 0.622 * 0.8 * VET (TEXP) _ Q ^ 2 2 * 0.622 *VPT(TCAB) 

PCAB - 0.8 VPT (TEXP) * VCAB - VPT(TCAB) 


( 6 ) 
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where 


TEXP = 86.9 + 0.066 * TCAB + 57-4 * (7) 

and VPT is the functional relationship between temperature and vapor pressure. 
VPT(TEXP) = 3207.0 * 10^ 

where fi = g** 60 - * * 3.244 + .009868 

* (647.27 - IS 5^2)) + (1.17 * lO" 8 ) * (647.27 - ) 3 (8 ) 


QCOHD(l) = 3.04 * (T(l) - T(2)) 

««««», -385 * .04l8 * PCAB * 144.0 * EM * CPG 

QBSEN1 = 48.3 * ( TCAB + 46o) 

*((0.385^(1) + 0.086 * T(2) + 0.0287 * T(3) + 0.238 * T(5^ 
+ 0.2615 * T(6)) - TCAB) * (1.0 - 0.33 * (14.7 - PCAB)) 


(9) 


( 10 ) 


The remaining core temperatures (nodes 5, 9, 13, 17, 21, 25, 29, 33, 37) are 
computed by similar equations. 

The muscle compartments introduce a set of nonlinearities in the blood flow 
calculations and metabolic heat generated with the shivering equations and 
blood flow regulations . In the example of the trunk muscle : 


QMET (6) = QB (6) + 0.3 * (UEfP * (PM-QB)) + O.85 * Q3HI7 (ll) 

where 

QS3HIV = .( -50881.0 + 1928.6 TC + 61-.4 TC 2 
+T(42) * (4453-4 - 168.5 TC + 5.37 TC 2 ) 

+T(42) 2 * (-106.97 + 4.05 TC - 1.89 TC 2 ))* 3.6 * 3.97 (12) . 

and TC is (T(l) - 32) * 1 

9 

QC0Nv(6) = bf(6 ) * (T(4i) - t(6 )) ■ (13) 

.There BF(6) = 13.2 * °-3 * (1-0-UBF « (BMgASAL)) + 0 .g 5 * ^ (lU) 

x • o 
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The elements of the skin have some additonal heat exchange terms as given 
in the heat balance for the skin of the trunk. 

dT(8) _ 1 *(Q^ET(8) + QC0HD(7) - QC0NV(8) - QRAD(8) - QSEir(8) 

dt C(8) 

- QhAT(8)) (15) 

where 

QRAD(i) The heat loss due to radiation distribution to the i^* 1 
element . 

QSEN( i ) The sensible heat loss due to evaporization of sweat in 
the i 6 element. 

QLAT( i) The latent heat loss due to evaporization of sweat on 

the i tn element. 

The skin elements also add to the complexity of the linearization process 
due to radiation terms based on Stephan's Law: 

« = A**"* (I ^ kln -I i n ) (16) 

where a is the emissivity of the skin with logic for including the garments; 
the equation for radiative heat exchange for the skin of the trunk becomes: 

QRAD(8) = HR * (TUG.(8) - TW) (17) 

where 

HR = AR(8) * EUG((TUG.(8) + 46o) 3 +(TUG(8) + 46o) 2 * (tw + 46o) 

+ (tug( 8) + 46o) *■ (TW + 46o) 2 + (tw + 46o) 3 ) (18) 

where , • 

TUG(I) = HR * TW + HC * TCAB * ML ^ + ^ (1 9 ) 

CLO 

HC is the convective heat losses from the skin. 

HC “ 0.026 * AC * (PCAB * VCAB) 0 ' 5 (20) 

or 

HC = 0.06 * AC * (PCAB 2 * G * ItUG(8) - TCAbI ) 0,25 (2l) 


whichever is greater* 
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The resulting first-order differential equation of degree 4 contains 
eight variables of state. Further difficulties are found when consid- 
ering the blood flow rates and the production and evaporization of sweat 
to these compartments. The final equations can be found in Reference 3 
and again constitute nonlinearities of the controlled system which are 
necessary to represent the thermoregulatory system. This set of 42 
equations for the controlled system includes nonlinearities which do 
not allow a valid Taylor Series expansion, since these equations cannot 
be expressed as a power series, oo M (see Reference 4 ), There- 

X C xjr (X-A) 

M=0 


M 


fore, the equations of the thermoregulatory model are not amenable to the 
linearization method's found in Miniquasi. 


In conclusion, it has been found that the methods used in MINIQUASI 
cannot be applied to the complex equations of the thermoregulatory model. 
If somehow these nonlinearities could be removed, the resulting equa- 
tions would, most likely, no longer accurately represent the physical 
processes of heat flow in the human body. 
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